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A deterministic optimal design problem for the heat 

equation 


Heiko Gimperlein*' Alden Waters* 


Abstract 

For the heat equation on a bounded subdomain O of R d , we investigate 
the optimal shape and location of the observation domain in observability 
inequalities. A new decomposition of L 2 (R d ) into heat packets allows us 
to remove the randomization procedure and assumptions on the geometry 
of f 1 in previous works. The explicit nature of the heat packets gives new 
information about the observability constant in the inverse problem. 


1 Introduction 

This article considers an optimal design problem for the heat equation: What 
is the optimal shape and location of a thermometer if we would like to recon¬ 
struct the heat distribution in a domain? We aim to introduce techniques from 
microlocal analysis, related to heat packet decompositions, in order to address 
a rigorous formulation of this question. The explicit nature of the heat packets 
sheds light on certain randomization assumptions and technical hypotheses in 
previous works. 

Optimal design problems for the placement of sensors have attracted much 
interest in analysis and computational mathematics. Recent works include 
[3 El 1 E H2, which consider observability and optimal design problems for 
wave and Schrodinger equations. For the heat equation, an observability esti¬ 
mate goes back to m, and P3] considers a simplified optimal design problem 
with random initial conditions. The dissipativity of the heat equation makes 
the reconstruction of high frequencies much less stable, one of the issues ad¬ 
dressed in this article. 
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For a precise statement of a model problem, we consider a solution u(t, x ) 
to the heat equation ut = A u in a smooth domain with homogeneous 

Dirichlet boundary conditions and arbitrary initial condition u(0, •) € C£°(fi), 
c.f Appendix Theorem 14.11 for well-posedness estimates. Given T > 0 and a 
measurable subset w c 11, we denote by Ct(oj) the best constant such that 


Ct(uj) / \u(T,x)\ 2 dx< 


fj 


u(t, x)\ 2 dt dx 


( 1 ) 


when it(0, •) € C£°(fi). Ct(w) gives an account for the well-posedness of the 
inverse problem of reconstructing u from measurements over [0, T] X oj. 


This article studies for which subdomains the heat equation is observable, 
and whether there are optimal ones. In the paper [9], a parametrix to the 
linear Schrodinger equation is built. We build solutions to the heat equation 
after their model, but we change the variable it e->■ t. We rescale their original 
high frequency data to obtain information about L 2 (0) data. Thus, a key new 
ingredient in our investigation is a decomposition of the initial data 


u = 




x + 2ij n t\ 2 \ 
4(cr 2 + t) / 


of u into heat packets following [9]. Here c n are constants, belongs to a 
lattice in M rf , and cr is a frequency parameter; we specify them later. While 
this decomposition is valid in whole space, we may reduce to this case after 
approximating the heat kernel on fi. The heat packets replace the propagated 
Dirichlet eigenfunctions of D, but still allow almost explicit calculations. 

We exhibit a family of optimal design problems which accurately approxi¬ 
mate the true optimal design problem when the initial data is not high frequency 
and the time scale is of the same size as the frequency. We let C^(uj) denote 
the constant associated to this approximate optimal design problem. We let 
0 < c <1 be a small constant; the precise conditions on e will be specified in 
terms of a fixed positive constant r/ with rj < 1 and the timescale T. The param¬ 
eter e describes the number of initial data points we need in our frame-based 
approximation. The parameter rj describes the convergence of the approximate 
observability constant C^(u) to C't(w). In particular we show that e needs to 
increase with the timescale T, and e may increase monotonously with r/. 

For computational clarity, we assume that i/j € with supp-0 C 

[—1, l] fi . Our initial data u(0,x) we say is of the form 



with eo a number in (0, ^diam(D)) which is sufficiently small enough for suppu(0, x) C 
0. The data is at minimal distance <5 from 50 with 5 € (-^/io, 1). We search for 
an optimal subdomain u in the set 


A4 m = {w C D | lo is measureable and of Lebesgue measure |w| = M |D|} . 
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It accounts for the fact that we measure the solutions on a subdomain of II 
with a fixed volume. 

The classical approach of }14] involves separation of variables using a basis of 
eigenfunctions ATj = —A j^j. Here one would decompose the solution into this 
basis as u(t,x ) = Yl < jLi a j e tXj ^j( x )- If we define bj = aje~ TX i , the question 
becomes to examine 
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dx dt 


: a 


/ g (\j+\k)T _ 


iTr f \ill 


where a denotes the spectrum of the matrix. This is a hard spectral problem 
since little is known about ’Pj(x)'k/ C (x) dx even in the case of the disk: the 
restriction of inner products of arbitrary Bessel functions to subsets uj C H 
cannot be computed explicitly. 

In order to avoid this problem, Privat, Trelat and Zuazua m replace a,j by a 
sequence of real-valued random variables and thereby introduce 

a random field u u (t, x). The /3J are independent identically distributed, of mean 
0, variance 1 and fast decay (e.g. Bernoulli). They then study the case of an 
averaged observability constant 


, ( e(A; — i f 

C r T and { w) = inf ° 1 A + (E/?,TO J dx 

/ e 2XjT - 1 f . 2 A 

= infer I --- / 'Pj(x) dx 

\ 2A j Ju 

JlXjT _ i r 

= inf ---- / T ; (x) 2 dx . 

j€N 2Aj L A J 

From the above we see that removing the randomization hypothesis reduces to 
a difficult spectral problem: To determine inner products on the subdomain uj 
of eigenfunctions for the Dirichlet Laplacian on H. We instead introduce an 
approximate problem, denoted by A which consists of a finite number of inner 
products of explicit functions. For the uniqueness results, we must randomize 
the initial data. We also assume that u(t,x ) intersects the observation domain 
uj at some time to £ (0, T) and this intersection is bounded below by mo- We 
set mi = mo||i/ ; |lc-o(fi) 7 b This criterion is related to [3j- 
Our main result is: 


Theorem 1.1. Let M € (0,1) andT < nhn{eg +s mi(C S! d||?/^(x)||c 2 + S (Q)) -1 , mid 4 }, 
s > d/2 an integer, and HcK^, then 

a) Ct(oj) > (uj) > 0 for all uj € Mm- 

b) There exists a unique uo* £ Mm such that C^’ rand (uj*) > C^’ rand (uj) for all 
uj € Mm- 

c) The subset uj* is open and semi-analytic. In particular, the Hausdorff mea¬ 
sure of duj* is 0. 
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The randomized initial condition required in pH] , which also included strong 
implicit assumptions on the geometry of 17. In order to study the randomized 
problem further, Privat, Trelat and Zuazua require strong assumptions on the 
level sets of these eigenfunctions ( (Hi) in pH]) and thereby on 17. Here, we show 
that removing randomization amounts to the computation of a finite number 
of inner products of explicit functions, for which the assumptions on level sets 
can be verified. We show that the problem can be modeled for arbitrary initial 
data in L 2 (!7) by Gaussians, suggesting that a deeper analysis of the heat kernel 
could be useful in including the randomized terms for the uniqueness results. 
In particular we prove one cannot ignore the non-randomized terms entirely 
in this scenario. We provide quantitative upper and lower bounds on their 
contribution to the observability constant. 

The number of grid points required for the reconstruction of C^(u) in¬ 
creases drastically at high frequencies. Of course, the initial data must be 
in iPg(17) for the solution to even exist, so in some sense this phenomenon is 
to be expected. The poor convergence for the high frequencies suggests that 
the use of monochromatic waves, which is introduced by randomization intro¬ 
duces a significant loss of information in the optimal design problem. In par¬ 
ticular, the randomized observability constant only provides an upper bound 
C^ and (oj) > Ct(w), while we have 

C't(w) > C$(oj). 

In essence, we are looking at the worst case scenarios for the behavior of ui. 
The reason we cannot extend to infinite times, is that the errors in the approx¬ 
imation argument dominate for large times. We also give an example of high 
frequency data showing the difficulty of removing the randomization assump¬ 
tion of [T21 . 

It should be noted that m is the first paper to construct an optimal observ¬ 
ability domain lj in any sense, and we rely on their ideas extensively. It is an 
earlier result of [10] that Ct(oj) > 0 for compact connected smooth manifolds. 

The heat packets lead to a numerical approximation scheme, where one 
determines a subset u* N € M.m with the best constant in ([U) among finite 
linear combinations of heat packets 

u — ^ ^ C n (f) n . 

\n\<N 

A basic problem is the convergence of lj^ to w* as 7/ —» oo in Hausdorff 
distance. There are well-known counterexamples due to spillover phenomena 
[8j for hyperbolic equations. However, the penalization of high frequencies in 
the heat equation allows to prove such a result, similar to fb4j for randomized 
initial conditions: 

Theorem 1.2. a) uj* n is uniquely determined, semi-analytic and open, 
b) The sequence u* N converges to oj* as N —>• oo. In fact, = ui* for large N. 
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Section [2] considers the auxiliary problem in the half-space, in which explicit 
computations based on the Gaussian heat packets are possible. The model for 
C^{oj) and the proof of Theorem 11.11 a) is given Section [21 The proofs of The¬ 
orems [TTT] b),c) and ll. 2 l are established together in Section [3j The construction 
of the frame implies that we have precision over N. For us, the results indicate 
N increases with an increase in T. 

Notation: For z € C d we write \\z\\ 2 = z T z for the square of its length and 
| z | 2 = z T z for the analytic extension of the absolute value on M. d . 


2 Truncation to model problem in upper half plane 


We consider our initial data to be a smooth compactly supported function, and 
explain how this can be used to obtain information on optimal design problems 
for more general L 2 (M. d ) data later in the paper. For general conditions, the 
evolution is approximated by a superposition of heat packets. We start with 
the upper half plane as our model and use the Feynman-Kac principle to move 
the problem to the bounded domain P. 

We assume diam(P) 77 , and |P| = 1. We extend 7 /' by 0 to a smooth 
function on M d and use the following lemma to decompose the function into 
heat packets. For the computation using Kac’s principle it is convenient if the 
origin, without loss of generality, has the maximal distance -ydiam(P) to <9P. 

We decompose ip as follows: 


Lemma 2.1. [Generalization of Lemma f.lf in j£|// Fix ip € C^°(M. d ) as above 
and 77 € (0,1) and eo € (0, ^diam(P)). We define ip eo {x) = e[[ d ^ 2 i[’(x/eo). Let 
<5 denote the distance of the center of ip eo (x) to the origin. Then there is a 
small 0 < e < 1 depending on 77 , eo such that the following holds: For a = 
y/eo5 log(e _1 ), L = uloglog(e~ 1 ) there are coefficients {c n } ngZ d and fn = x 
with: 


I / \ c ( \ x \ . j . 

Mx) - L exp r 4^ + 


x 


<V\\lpe 0 {x)\\ L 2 {u) . (2) 


L 2 H 


Moreover, 


\c n \ <^,k C k ^(ae 0 L 2 ) d/2 min jl, (e 0 |£ n |) fc | 
for all k S N. One may choose c n = 0 unless 

n g 5 = {t 6 Z-i : „ 1 o 8 4 (< - 1 , < j&| < 
if we include all the vectors then rj = 0 . 

The proof in [9] is for d = 3. We include it here to show that the smallness 
condition on e is compatible with later restrictions on e in our paper. The 
proof in [ 2 j is for high frequency data only, which we rescaled to obtain our 
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result, since eo is arbitrary. The rescaling affects a and L. This implies the 
set of required vectors S becomes very large as eo —>• 0. The integer k reflects 
the regularity of if, and it is possible to bound k by d provided e is small, 
and still produce an (? sequence c n . The construction of the £ 2 bound requires 
eo < <5 < 1, but it should still be possible to construct an i 2 sequence as long 
as b > 0. 

Proof. We define 

c n = (27 xa)~ d J ?/> eo (x)(27rcr 2 ) d/4 exp - in ■ dx 

[—7rL,7rL] d 


and 

4>n(x) = (27T(T 2 )' d/4 exp ^-^2 + * n ' ' 

Note that 11(3?) 11 z, 2 (M d ) = 1- As the Fourier series of a smooth function con¬ 
verges uniformly, we have by Plancherel’s theorem: 

Ao(x) = c n (j) n {x) (3) 

nez d 


on [—7rL,7rL] rf . If e is sufficiently small, supp^ eo C [— The series (0) 

gives the result immediately if uo C fl, which it is. The error comes from the 
truncation, which we will prove a bound on using the estimates on the \c n \. 

We now prove an upper bound on the coefficients. Let M\ be a constant 
such that ||^||c d (r2) — Tfi. From the definition of c n we obtain 

~.d/2 /1 i2 \ ^ 

\c n \ < -^||^e 0 (^)exp (^J|| L i( R «i) < • 

In order to derive the upper bound, we use integration by parts. Let B = 
i yhp • V. We see that B fc exp (—in ■ = exp (—in ■ fj. The adjoint of B is given 

by B* = —iV • |hp-. We therefore obtain the following 


\c n \ = (27 tL) 


= (2irL) 


-d 


—d 


J B fc exp (in ■ if eo (x)(2T\ rcr 2 )^ 4 


exp ( ^-2 I dx 


s 


exp i —in ■ ^ 


L 


-o d/2 ' t l ) ta{x)(2'KG 2 ) d/i exp (^ 2 ) 


dx 


< L 


—d 


T \ k / \ 

L \ f a \ 


\n\ 


d/2 


e oJ 


E I' 8 " 


|a|</c 


* 1 J ) exp (§0 


^k 


M d / 2 ( L 
L d y €q |n 


Ml 
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The claimed upper bound follows. For the remaining claim we notice that 



Fix N 6 N. If n < N, we use the upper bound for c n to estimate 

_2 


^ ^ Cn4 > n(.'E 

|n|<7V 


(ceo) d \ ^ 

~ -JIT- ex P 


2L? 


\n — m\ 


< 

r^j 


(¥)' 


Mi 


|n|,|m|<JV 

For n > N, we use the second upper bound for c n with k = d to obtain 


^ ^ Cn^r, 


|n|>JV 


d / r \ 2d 


< 


(o-e 0 ) / + 


L 2 (R d ) 




£ 

|n|>|m|>Af 


1 1 

n| 3 Iml 3 


exp 


<j- 


2 L 2 


n — m 


< 


~) 


) d E o 


< 


|n|>|m|>7V 


L V 


2d ~ V eoiV,) 


Mi . 


We conclude 

^ ^ Cn4>r 

w<- 


+ 11 CnAill^QRd) < (loglog(e *)) d M 1 


M>^l°gl°g(e X ) 


: £Q lo S lo g( E H 

The assertion follows provided that 

(loglog(e^ 1 )) _d Mi < r] . 


(4) 


□ 

This observation of [9j gives a precise expansion of the initial condition into 
Gaussian heat packets 

<W*) = {27T ^ )d /A eX P (- |X 4a 2 01 + - X °)) ’ 

which evolve in according to simple analytic formulas and are centered 
around a suitable xq. In the decomposition we selected, xo = 0, by our choice 
of origin. Here xo denotes the fact that we could have had more that one ip eo 
in our decomposition for u(0,x). We leave xq in to show that it is difficult to 
do the Gramian computations when the centers are not suitably separated. 

For simplicity we write 


4>n,x 0 (t,x) = [ex.p(tA R d)(j) n , X0 }(t,x) , 


where we have that 


&n,x o (t,x) — ( 


a 


\V2tt (a 2 + 1) 


d/2 


exp 


\x — x 0 + 2i£ n f| 2 
4(cr 2 + 1) 


^|£n| + 'i’Cn ' {x Xq) 
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We let 


k Rd (t,x,y) = (Ant) d ^ 2 exp (—\x — y\ 2 /At) 


denote the heat kernel in M. d and kF the corresponding heat kernel in fi. We 
take Kac’s principle, Proposition 6.3.1 in |2], to approximate the evolution of 
the heat packets in the upper half plane. Kac’s principle states: 

Proposition 2.2. Assume that Pi C. M. d is smooth, and let x € fi he fixed. For 
every y G we let 


It follows that 


to(y) 


d(y,T) 2 
2 d 


0 < k Rd (t,x,y ) - kn(t,x,y) < < 


(Ant) 2 exp(— d(y, T) 2 /At) if 
(Ant 0 (y))~^ exp(-d/2) if 


t < t 0 (y), 
t > t 0 (y). 

(5) 


A short computation using Kac’s principle shows that: 
Lemma 2.3. We have for 0 < t < y o<5 4 , and HV’Hc'Offi) < Afy. 


llexp^An)^^)-exp(tA R d)^ eo (a:)|| c -([ 0jT ] X n) < yM 0 

whenever e$ < 5 2 < A. 

Proof. We would like an estimate on 


sup | / (k Rd (t,x,y) - kn(t,x,y))fj eo (y)dy \ . (6) 

t,x J 

R d 


We divide Pi into a region I, defined by d(y,dPl) > t 1 / 4 , and its complement, 
region II. Using Kac’s principle we can bound ([6]) from above by 

sup J (Ant)~ d / 2 exp ^ \if eo (y)\dy 

I 

/ [ 27T \ 

\~d) ( d (y> dQ ))~ dex P(~ d / 2 Me 0 (y)\dy . 

ii 

We introduce coordinates near the boundary dPl with y\ = d(y, dPl) and y' = 
(j/ 2 ) • • •, yd)- For the integral in region I, we have after change of variables 


diam(f2)/eo 


f d 

e 0 


(Ant) 


-d/2 


exp 


_ 7 , 2,2 

Vi 6 o 

At 


\F(yi,y')\dy'dyi < 


d exp 


i 

'ip72 


: 0 


(Ant) 2 


M 0 \dPl\ 


i 1/4 /eo V'/ e o 








In region II the integral can be estimated by 


£1/4 


2-7r\ _ f 

-j-j exp(-d/2)tq d J, |</k 0 (di, y')\dy' dyi 


i (2ir\~ d/2 


= £n 


d J 


exp(—d/2) 


y 

t 1/4 A o 

/ 

5/eo 


Vi d \^{yi,y')\dy' dyi 
Jy'/to 


=4^v^xpm / 2 ) m » £ °"- 


d / 


d- 1 


-\dn\ . 


and the integral is empty if t 1 / 4 < 5. If we impose this condition, then quick 
inspection of the second integral gives the desired result. □ 


Initial data closer to the boundary would result in the construction of 
a reflected heat packet, which we leave as an area of improvement in our 
analysis. For the construction of the parametrix, we impose the condition 
|w|Mor/o < rj\\u(to,x)\\ L 2 r UJ ) for some to € (0,T). Because we are dividing by 
the measure of the set |w|, for r/o to exist, we impose the condition similar to 

m : 

(•) The support of the solution u(t, x) must intersect the observation set oj at 
a point xq for at least one time to, and |ri(to,£o)| > mo- 


The condition implies rjo < as m(w)Mor?o < t i\\ u {Q, x )\\l 2 (uj)- 

Lemma 2.4 (Analogue of Proposition 4.7 in [9]). The evolution near the bound¬ 
ary is given by 

exp(tA Q )?/> eo (x) = 5><^ 0 (t,x) + v(t,x) , 

n 

where 


\\v(t,x)\\ c o {[0 )T] . L 2 H) < v\\M x )\\l^( u ) ■ 

Proof. This estimate requires two steps. For the first step we have: 

II exp(fA Rd )£/ eo -exp(tAn)^ eo ||z,2 ((j) < ^||V'eo(*)IIz, a («) 

by Kac’s principle and (•). For the second step by inequality (ED and the 
parabolic maximal principle 


L 2 M 


exp(tA R d)V’ eo -E 

n 

exp(tA K d) ^ eo -E CnfirijX o 


l 2 h 


— || exp(tAjjd) \\l 2 (uj)^-L 2 (uj) HV’eo ^ ] C n (f> n ^ xo (x) ||l 2 (ij) 


< |ll^eo(®)llL 2 M , 
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c.f. Q], Theorem 1.3.3. 


□ 


In particular, we notice that away from this regime we can approximate 
the evolution of the heat packets by the evolution in M. d . This approximation 
simplifies the analysis considerably, so we choose to focus on this regime as our 
model case. 

The key ingredient will be precise estimates for the inner products of heat 
packets on M rf , which follow from the explicit Gaussian shape of <f> n ,x 0 - Analo¬ 
gously to m , we consider the Gramian matrix G corresponding to the evolved 
heat packets: 


G nm (x 0) HO) ^0 



(fi •^')4 > m,yo (t) x) dx dt . 


Inspired by techniques used for Gaussian frames in [18] . we would like to show 
that the largest contribution to the observability constant comes from the di¬ 
agonal terms, and that the largest contributions are due to low-frequency heat 
packets. In m a randomization assumption was required to remove the off- 
diagonal terms. 

To simplify the notation, set 


Aft, n ) = exp 


2hx 2 |£ n | 2 \ 

a*+t J ' 


We define B(x o, V + t) as 

B(x o, a/ct 2 + t) = {x : \x - x 0 \ < 2(t + cr 2 ) 1 / 2 } . (7) 

Proposition 2.5. Let e\ < 1 and T = ei<r 2 . We assume e in Lemma, \2.1\ is 
small enough such that 12 C [— a 2 , a 2 ] d . We let Cd be a fixed constant depending 
on the dimension only. We have the following estimates: 


a) Lower bound for the diagonal terms: 

G nn {x o,x 0 ,uj) > C d (exp(-l))a d f ^ ^== A(t,n) dt . (8) 

Jo \B{X 0 , \Ja A + t)\ 

b) Upper bound for the diagonal terms: 

G nn {xo,xo,u) < C d a d (l + erfc(l))(erf(l)) _1 [ + ^ A(t,n) dt 

Jo \B{x 0 , \Jo A + f)| 

(9) 

c) An upper bound for the off-diagonal terms: 

c mJ o A 1/2 (t, n)A 1/2 (t, rn) exp ( - ^ ) dt (10) 
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d) An equality for the off-diagonal terms when |n|, \m\ < L, a 2 + t < l, u is 
radially symmetric with suppw C [— R, R] d : 


| G nm (X(), X(), W) | 


(11) 


= Cd\aj\eifR 


/ o 


cr 


a 2 +t 


d/2 


1 + (o(l)(a 2 + t)y) (A 1 ' 2 (t, n)A}G(t , m)) dt . 




Proof. The estimates follow from the explicit formula: 
/ / <t>n,x 0 {t,x)<j> m!X0 (t,x) dx dt = 

Jo Ju d 


( 12 ) 


exp 


2(er 2 + 1) 


T *X • I (fn f,Tn) 


a 2 + t 


exp -(£n + £m) 


ta l 


a 2 +t 


Jl 

Jo J K c 

In a) we integrate the Gaussian over the intersection of uj over the standard- 
deviation part of the integrand around xo- We recall the following definitions 
of the error and complementary error functions: 

oo b 

/ exp(—x 2 ) dx = erfc(6) , —= / exp(—x 2 ) dx = erf(6). 

V7T J 

b o 

The error function has the following asymptotic behavior [1]: 

1 


erfc(6) 


Vb > 0 


(l + ai + a 2 b 2 + a 3 b 3 + a 4 6 4 ) 4 

with ai = 0.278393,02 = 0.230389,03 = 0.000972, a 4 = 0.078108. 
For part a), we start by integrating 

f-Oin(xO) Xq, R ) 


(13) 


in x, and notice that 


B 


exp 


V a 2 ~et 


Vcr 2 +t 


~ erf(l) . 


Then we may replace the limits of integration in the inner product over 
to those of the ball B(x o, V& 2 + t), with small error as dictated by (fT3l) . The 
desired result follows by a change of variables in x and using 

infexp(—x 2 ) = exp(—1); 

B 

here the limit is over the rescaled ball. The result in part b) follows similarly. 
For the result in part c) we see 

d/2 

X • • • 


/ / (j>n,x 0 (t,x)(j)m,yo(t,x) dx dt < 

Jo J R d Jo 

i2 (l^n| 2 + |£m| 2 ) (x 0 — 2/o) 2 


(T- 


x exp 


a' 2 +t 


8(<r 2 + 1) 


a 2 + t 

- K(i n - e 


- t(\Zn\ 2 + \U\ 2 ) 


dt 


(14) 


dx dt 
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We notice that exp(— K{^ n — £ m ) 2 ) is decaying for short times. Unfortunately, 
we cannot easily exploit the decay in K. In particular, when integrating with 
respect to x in (fT2D . this amounts to an estimate on 


J exp(ifc • x) exp(—ax 2 ) dx 

UJ 


(15) 


for some k, a constants without losing information on the oscillatory part. By 
Hardy’s uncertainly principle, c. f. na Theorem 1, the only way a Gaussian will 
Fourier transform to another Gaussian is if the initial data is Gaussian, and we 
are integrating over all of M d . Otherwise, we obtain a larger rate of decay. We 
define a smooth cut-off function Xu, which is 1 on the support of u>. We see 
that 


Kx^,exp(-ax 2 )exp(ifex)) i 2 (K d ) | < 

IIXw|li 2 (s 0 )|| exp(—ax 2 ) exp(zfcx)|| 2 2(R d) < C d m(uj)a~ d/2 

with C independent of a and k, while 

r2 k 2 ' 


. d 
7T\ 2 

- exp 
a / 


7r 


= / exp(— ax 2 ) exp(ifcx) dx. 


Thus, we can bound |G nm (xo, xq, w)| with a loss of exponential decay. It follows 
that ([HD is bounded by 


C d m(uj) 


rT / a 2 \ d/2 


ex P I — (Cn + Cm)' 


t(T~ 


, 0 \<? 2 + t) ^ v vsn 1 sm v 2 + ty exp v 4(u 2 + t) 

Recall as a modification of Lemma 3 in m 

Lemma 2.6. Suppose |/(x)| < exp (—ax 2 ) for all x in {x : |x| < R} and some 
a > 0. Then /(£) is smooth and one has the bound for £ £ / : 


(To - yo? 


dt . 


j£P 

'a \ a 


for any interval I, such that |/| < 1. 

Proof. The Fourier transform has the Taylor series expansion 

/(C) = /(0) + £<%/(°) + y<%/(o) + Odd 3 ). 

We can then calculate 


^/(C)l$=o = J (27rix) fc exp(— ax 2 ) dx 
\x\<R 


(16) 


(17) 


(18) 


□ 
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Our function which is the integrand in (1121) satisfies the criterion of Lemma 
in Id. We apply the Lemma which generalizes easily to W l . It follows from 
the equality (1111) the off-diagonal terms contribute substantially to the matrix 
norm of a Gramian if cr 2 < 1 and |n|, \m\ < L. Examining the leading order 
term part of the contribution in dill) , and the upper bound in c) we see 


erLRinf((A 1/2 (f, n)A 1/2 (t, m)))\\f) ntXQ (t, x)\\] 


C d 


l 


lL 2 ((0,T)xK4) — 


<7 


a 2 +1 


d/2 


, m )) dt < 


sup((A 1/2 (f, n)A 1/2 (t, m)))\\(/)n, xo (t, s)|li 2 ((0iT)xR d ) 

Both prefactors mf((A~ 1 / 2 (t,n)A 1 / 2 (t,m))) and sup((A” 1//2 (t, n)A l / 2 {t, m))) 

* t 

can contribute substantially when summing over n m. □ 

We need the following short time estimate. 

Lemma 2.7. Assume that the C 2+s norm of ip is bounded by M 2 with s > d/2 
an integer. Then for all 0 < t < o6q +s we have 

I u(t,x) -i/) eo (x)\ < r/o. 

with C Sj d a constant depending on d and s only. 

Proof. The mean value theorem states that 3to € (0,1) such that 

\u(t,x) - Ao( x )\ < \d t u{t 0 ,x)\t . (19) 

We know by definition dtu = A u, and by Sobolev embedding with s > d/2 

Il^' u llc°(f2) — C S: d\\u(t, x) ||i 72 +s(Q) < C Sjd \\u(fi, 3 ')II_H' 2 + S (n) 

The last line follows from standard energy estimates (Appendix, Theorem 14. ID 
and the finite expansion from Lemma [2.11 Our choice of timescale clearly works 
by scaling. □ 

We conclude from the precise values of the upper and lower bounds in 
Proposition 12.51 that: 

Corollary 2.8. Let T < e and e sufficiently small. We have the following 
bound 


V Jo InXu(x)\u(t,x)\ 2 dx dt 
2 \u(T, x)\ 2 dx 




n'ffin^x 0 


(t,aP" 2 


< 


L 2 (wx(0,T)) 


< 


S Cn^n^x oC^\ x ) II l 2 (Q) 

n v ' 

2 fo T fnXuj(x)\u(t,x)\ 2 dx dt 
g f n \u(T,x)\ 2 dx 


( 20 ) 


with C^(oj) < Ct{oj), and C^(lo) > 0. The replacement is possible even if we 
consider a finite number N of (j> ntX0 (t,x) with |n| < N as dictated by Lemma 

lOl 
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Proof of Theorem U. 1 \ a) and Corrollarv \2.8\ . . We recall that ||<^n.o;o(Oj ;c )llf, 2 (M d ) = 

1. Recall that 77 is the small frame parameter given to us by (|2|) . We make the 
natural assumption which we call the bootstrap assumption, 

K t,x) - u(0, x)\ < 770 =>- ||u(t,x)|| Z/ 2 (w) > 2 i]\\u( 0 ,x)\\ L 2 ( U} ) , 

which is valid for short times and small 770 such that \u\r]Q < rj\\u(to, x)\\ L 2 ^ 
for some to £ (0, T). Because we are dividing by the measure of the set |w|, 
for t/o to exist, we must have the criterion (•) hold. The question of validity of 
this assumption for arbitrary eo is given by Lemma which fails as eo —>• 0 
because the timescale shrinks with eo- We also know 

\W(t,x)\\ L 2 ^ < || exp(fAn)(u(0,x) - E Cn 0 (0, ^))||l 2 (u;) + II ^ ^ C n (f>n,x 0(^5 ^)||l 2 (o;) 

n n 

< 77|| 17.(0, x)|| £ 2( w ) + \\^2c n (j) n:XO (t,x)\\ L 2 {uj) . 

n 

By the bootstrap assumption and the energy estimates in the Appendix, we 
have 


IWM)I|l 2 (u,) -» 7 |KO,x)|| £ 2 (w) > ??||u( 0 ,x)|| i 2 H > 77||7i(t,x)|| £ 2 (w) . 

We see also that 

|l 5 ^c r ,^ r , iXO (r,x)|| £ 2 ( n ) < \\u(T, x)|| £ 2 (fj) + r]\\u(0, x)|| £ 2 ( q) < 2||u(T, x)|| £ 2 (n) , 

n 

again by the bootstrap assumption and the energy estimates. The desired result 
(|20|) follows. We know from pjQ] that the left hand side of (|20l) is nonzero, so 
C^(u) is nonzero. 

Therefore we have made only the following assumptions on e in terms of the 
parameter 77 : 

0 < (loglog(e“ 1 )) _d Mi <C 77 . 

This was used in Lemma 12.11 on the condition (|2J) precisely dictated by dU) . We 
also made the short time assumption in Lemma 12.71 and condition (•): 

t < Cf}r]oel +s M 2 -1 , M 7707770 < M 0 r]\\u(t 0 , x)\\ L 2 ^ for some t Q € (0,T) 

and the assumption 12 C [—cr 2 ,(T 2 ] d for the bounds a,b) in Prop [2751 to hold. 
We notice that the short timescale is only necessary for the bounds on little 
u). For longer timescales, the parametrix itself is valid on 12. The timescale 
for the parametrix to be valid is largely dictated by Kac’s principle and can be 
improved by including reflections at <912. □ 

Remark 2.9. For an arbitrary collection of data that there is no interaction, 
e.g. G nn (xo,yo) can be not that much different from G nn (x o,xq) unless there 
are some strong hypotheses on the separation of xo,yo- Also, the frequency 
scales of the relative Gramians change as L and a scale with the respective 
distance to the origin. We leave this as an area for improvement in our analysis. 
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3 Proof of Theorems 11.11 and 11.2 


We re-label where it is understood so that we are now studying the randomized 
initial field, e.g. c n i->- /3^c n . We drop the subscript xq where it is under¬ 
stood. The estimates on the randomized field are for arbitrary bounded times 
T < a 2 , and we only need the short time assumption for the approximation to 
randomized observability constant. 

From the previous section, we conclude that our appoximate randomized 
obervability constant with heat packets can replace the observability constant 
with u u a randomized field as in the Introduction. However, we cannot obtain an 
appoximation of the true observability constant because the off diagonal terms 
contribute substantially to the matrices required, by Proposition 12.51 c.d). In a 
deterministic problem, because the heat packets are not orthogonal, one obtains 
a generalized eigenvalue problem Gc = XHc, which only becomes an honest 
eigenvalue problem if one discards that H is not diagonal, i.e. if one neglects 
the L 2 -inner products between different and assumes they are orthogonal: 
H nm = f n 4> n (T,x)4> m (T,x) dx ~ C n (T)5 nm . All the analysis below is only 
possible if H nm = f n cj) n (T, x)cj> m (T, x) dx ~ C n (T)5 nm , whence the need for 
randomization. Without this assumption, one is forced to deal with Gc = XHc, 
and an analysis of both G and H is required: The relevant matrix is then G 
in the basis eigenvectors of H, and one needs that this is dominated by the 
diagonal. In the previous section, we were able to reduce the question to a 
finite matrix optimiztation with a lower bound on the observability constant. 
One could obtain an poor upper bound in terms of only the diagonal entries, 
but this bound is not very sharp. 

We need the following proposition to aid in the computations. Let Um = 
{Xu ■ w € -Mm} the set of characteristic functions supported on sets in A4 M . 
We study the optimal observability constant as a functional on Ul- 


Ct{Xlj) = inf 


Jo f n Xu>(x)\u(t,x)\ 2 dx dt 
f n | u(T, x)\ 2 dx 


( 21 ) 


where the infimum extends over all solutions u € of the initial-boundary 

problem for the heat equation with homogeneous Dirichlet boundary conditions 
and initial condition u(0, •) € C£°(fl). 

To assure existence of minimizers, it will be useful to study a relaxed prob¬ 
lem, in which we extend Ct(Xu) from Um to its closure in L°° with respect to 
the weak* topology, 


U M = {a£ L°°(H; [0,1]) : J a(x) dx = M\Xt\} . 

n 


We set: 


Ct{cl) = inf 


lo f n a{x)\u{t,x)\ 2 dx dt 
f n | u(T, x)\ 2 dx 


( 22 ) 


In order to understand the existence and properties of which maximize 
Ct, we would like to replace u( 0, •) by a superposition J2 n c nfin of heat packets 
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as in Lemma EU Then we define Ct(Xu) as the infimum over all admissible 
choices of c n . 

From the discussion above we focus on the analysis in the upper half plane. 
As in |14j . we may write this as an eigenvalue problem, i.e. with constraint 

= i. 

n 

However, here the time-dependence is not entirely trivial. From Corollary 
I2.8h l and b) we may then conclude that Ct{Xw) > 0, he. P ar t a ) °f Theorem 
11.11 Our proof of Theorem 11.11 and Theorem 1.2 will now follow closely along 
the lines of [14] . The basic hypotheses on the spectral decompositions there 
can be proven for our heat packets, but we retain the ordering (H 1 ) and {H‘>) 
from M- These seemingly natural hypotheses imply strong assumptions on 
the level sets of the eigenfunctions for the Dirichlet Laplacian which are never 
actually proved in [IT . 

The following Lemma will be used to show that minimizers of the relaxed 
problem are characteristic functions: 

Lemma 3.1 (Hi). Assume that there exist a subset E C 11 with \E\ > 0, an 
integer N £ N*, coefficients otj € M+, \j\ < N and C > 0 such that 

T 

otj / \<(>j(t, x )\ 2 dt = C 
\j\<N 0 

a.e. on E. Then C = 0 and otj = 0 for all j. 

Proof. The functions 

wt ' x)? = (tSFTT)) exp { JX ~ X S~tT V ) exp( - 2|f/ ‘) 

extend from to holomorphic functions of x on C d . Integrating in t , also 
T 

YO otj f \<fj(t, x)\ 2 dt admits a holomorphic extension and is constant on E. 
\j\<N 0 

Because \E\ > 0, E contains an accumulation point in Cl, and therefore 

T 

yy otj / \<f>j(t,x )\ 2 dt = C 
\j\<N 0 

is constant for all x € C d . We integrate both sides of this identity over x € 
to conclude that 



2ta%\*\ 

o 2 + t J 


dt 
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is infinite, whenever C > 0, a contradiction. Therefore C = 0, and therefore 



a 2 + t J 


dt = 0 . 


As all summands are nonnegative, and the t-integrals positive, we conclude 
(ij = 0 for all j. □ 


Let 

dj = ( j \(j)j(x,T)\ 2 dx 
\n 

In order to assure the existence of a solution, we use the relaxation as defined in 
g]. Because the set Um is not weak-* compact, we consider the convex closure 
of Um in the weak-* topology of L°°, which is then 

Um = {a € L°°( 0 ; [ 0 , 1 ])| J a{x ) dx = M\Q\}. 

n 



This relaxation was used in un nn ns]. If we replace Xu € Um with a € Um, 
we define a relaxed formulation of the optimal shape design problem, by 


sup J(a), 

O-Gt/ M 


where the functional J naturally extends to Um as 


J(o) = inf dj j a(x)|(/>j(x )| 2 dx. 


We show the following: 

Lemma 3.2 ( H 2 ). For all a € Um one has 

liminf f f a(x)d n \(j) n (t, x)\ 2 dx dt > 7 i(T) . (23) 

H-»°° Jo Jo. 

Here 71 (T) is the “appropriately renormalized,” first-frequency heat packet 

m Jo L\Mt,x)\ 2 dxdt 

71 / \(fi{T,x)\ 2 dx 

n 

Proof. It suffices to prove the estimate for Xw £ because a is in the weak-* 
closure of Um- The estimate (1231) is equivalent to 

.CL \fi(t,x)\ 2 dx dt Jq f w \cj) n (t,x)\ 2 dxdt 

fn |<()i(T,x)| 2 dx |n|— >00 f \f n (T,x)\ 2 dx 

n 
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We use the shorthand B t for the ball defined by (0). Let mf = and 

Cb = exp(—l)(erf(l)) as in Proposition 12.51 a.bh The left hand side has the 
upper bound 


{rr%)- Cd- 


(sup mf) f A(t,l)dt 

o _ 

C b A(T, 1) 


n\-l^ , ex P( 2 l6| 2 ^) - 1 


< (m“) C d C B VI 


2iai J 


< oo 
(24) 


and the right hand side has the lower bound 


(inf mt) f A(t, n ) dt 


o mpr'CdCs- 


> (m^)~ 1 C d C B 


A(T,n) 


M exp(2|£ n | 2 T) — exp(3/2|£ n | 2 T) 


\B t \ 


3/21 ^ 


(25) 


where we used the fact 

UJ II • p LJ 1^1 

supm ( = I a; | mi m t = — — 

t 1 

The right hand side of (|25|) goes to oo as n —> oo and the claim is verified. This 
computation can also be used to show C^{uj) > 0 directly, provided we recall 
the denominator in the randomized observability constant is bounded below by 
Y \ c n\ 2 by energy estimates in the Appendix. □ 


Standard variational arguments assure the existence of a unique relaxed 
solution: 


Lemma 3.3. The optimal design problem admits a unique solution a* € Um- 

Proof. As a result of Corollary 12.81 it suffices to examine the minimization prob¬ 
lem for the functional J on Um , 


J{a) 


inf 

cg£ 2 (N d ) 


Jq J n a(x)\ Ycn4>n(t,x)\ 2 dx dt 
n 

In Y I c n (j)n(T,x)\ 2 dx 


If we consider the normalization f Q Y I c n 0 n (T, x)\ 2 dx = 1, then it follows that 


we may equivalently consider 

rT 


inf / / a(x) d n \<j) n (t, x)\ 2 dx dt . 

Jo Jn „ 


For a € Um , the mapping 


a i->- d n 



a(x)\<j> n (t, x )\ 2 dx dt 


is linear and continuous in the weak-* topology of L°°. Therefore J is upper 
semicontinous as the infimum of continuous linear functionals. Because Um is 
compact in the weak-* topology, the result follows. □ 
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Proof of Theorem 1.1 b), c). We now define the truncated functional 

Jjv(o) = inf d n / a(x)\<t> n (t, x)\ 2 dx dt 

M<w Jo Jn 

for all a £ Um and consider the problem 


sup Jn (a) . 

Ul 

The above problem has at least one solution a N € Ul, by the arguments 
that proved Lemma 13.31 because Ul is weak-* compact and Jn upper semi 
continuous. Using (Hi) and (H 2 ), we are going to show that the solution a N is 
actually a characteristic function of a set c v N , a A = Xu N € Um■ Namely, if we 
denote 


&N — — (o ; n)|n|<Af • ^ ^ &n — 1} , 

\n\<N 

the Sion minimax theorem implies 


sup min d n / a(x)\6 n (t, x)\ 2 dx dt 

a eU M 1 - n - N Jo Jn 

= max min / a(x)tpN(x) dx 

a&U M Jn 

= min max / a(x)^N(x) dx . 
aeS N a <dU M Jn 


Here we have defined 



|<^> n (t,x)| 2 dt . 


As a result, there exists a N € Sn such that (a N ,a N ) is a saddle point of the 
functional. This then implies that a N is a solution to the problem 


max 
clGLJ m 


j a(x)tpN(x) dx . 
n 


By (H\), the functional (pN cannot be constant on a subset of H of positive 
measure. This implies the existence of X N such that a N (x) = 1 when <pn(x) > 
A n and a N (x) = 0 otherwise. As Jn is concave, the set of maximizers is 
convex. Since every maximizer is a characteristic function, a N £ Um must be 
the unique maximizer. We note that tpN is analytic and therefore ojn is open 
and semi-analytic. □ 

Proof of Theorem 1.2. It remains to compare the maximizers a N £ Um of Jn 
with the maximizer a* £ Um of J from Lemma 13.31 First, 


JN 0 (a*) < 7i (T) , 
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and ( H 2 ) applied to a* shows that 


T 

inf d n / a*(x) / \4> n (t, x)\ 2 dx dt > 7 i(T) (26) 

\n\>N 0 J J 

n 0 

for some Nq € N. From ()26l) we have that 

J(a*) = min / Jjv 0 (a*), inf dn [ [ a*(x)\<f> n (t,x)\ 2 dxdtX = Jn 0 (cl*) . 

{ M>iVo J 0 J 

If a N ° € Um is the maximizer of Jjy 0 , we show that J(a*) = JN 0 (a N °) as in [I HI : 
Indeed, because a N ° maximizes Jw 0 over Um, w e have that J(a*) = Jn 0 {cl*) < 
J No (a N °). We assume 

Jn 0 (cl*) < J No (cl N o ) (27) 

and obtain a contradiction: As Jjv 0 is concave, for every t £ (0,1] the assump¬ 
tion ([271) implies 

JN 0 (a* +t(a N ° - a*)) > (1 -t)J No (a*) +tJ No (a N °) > Jn 0 (o,*) = J(a*) . 

By the choice of No, one concludes that 

J No (a* + t(a N ° - a*)) = J{a* + t{a N ° - a*)) > J(a*) , 

in contradiction to a* being a maximizer of J. So indeed, Jn 0 (cl*) = J(a*) = 
J No {a No ), ot a* = a N °. □ 

Remark 3.4. Even though we are finding our optimal set over all indices n, we 
remark that the truncation of the admissible indices in Lemma 12.11 implies a 
bound on N For the case of high frequency data where eo —>• 0, IV increases 
enormously. 

4 Appendix: Well-posedness Estimates for the Heat 
Equation 

We consider the problem: 


dtu = A u 
u{t,x)\ xe9n = 0 
u(0,x) = g(x) 


(28) 


with g{x) £ Cc(Q). We claim: 

Theorem 4.1. Let A: £ No- Problem (1281) admits the following well-posedness 
estimate: 


u \\L°°(0,T);H k (tt)) + IMIl 2 (0,T);ff fe + 1 (fi)) — C\\g\\ H k(ty . (29) 
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Proof. We multiply (|28l) by u and integrate over Ll while applying the divergence 
theorem. We note that the boundary condition gives us: 


Id 
2 dt 


u 2 (t , x) dx + 


n 


n 


\Du(t, x)\ 2 dx = 0 . 


Integrating this equation with respect to time and using the initial condition 
we obtain 


1 

2 



t 



0 ft 


| Du(s, x)\ 2 dx ds 


1 

2 



dx . 


Taking the supremum over T gives the desired result for k = 0. The result for 
k > 0 follows by differentiating the heat equation and choosing \/ k u as a test 
function. □ 
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